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Proposal for exotic-hadron search by fragmentation functions 
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It is proposed that fragmentation functions should be used to identify exotic hadrons. As an 
example, fragmentation functions of the scalar meson /o(980) are investigated. It is pointed out 
that the second moments and functional forms of the u- and s-quark fragmentation functions can 
distinguish the tetraquark structure from qq. By the global analysis of /o(980) production data in 
electron-positron annihilation, its fragmentation functions and their uncertainties are determined. 
It is found that the current available data are not sufficient to determine its internal structure, while 
precise data in future should be able to identify exotic quark configurations. 

PACS numbers: 12.39.Mk, 13.87.Fh, 13. 66. Be 
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In the hadron mass region below 1 GeV, there are 
scalar mesons, /o(600), /o(980), and do (980), whose in- 
ternal configurations are not obvious [l[. Their fla- 
vor compositions could be /o(600) = (uu + dd) /y/2, 
/o(980) = ss, ao(980) = ud, (uu — dd)/\/2, ud in a 
simple quark model by considering the mass relation, 
m u ~ rrid < m s . However, this assignment implies a 
mass sequence, m(/o(600)) ~ m(a (980)) < m(/ (980)), 
which contradicts with the observed one, m(/ (600)) < 
?7i(a (980)) - m(/ (980)). If / (980) and a (980) are 
exotic states such as tetraquark ones, the observed spec- 
trum could be naturally understood. Since /o(980) and 
ao (980) are experimentally established resonances, they 
provide a good opportunity to study exotic mesons be- 
yond a naive gg-type quark model. 

First, a brief outline of recent studies is given for the 
/o(980) structure. In a simple quark model, a light scalar 
meson /o with J PC = ++ is identified as a 3 Pq quarko- 
nium with the flavor structure (uu+dd)/y2. However, if 
such an ordinary qq configuration is assigned for /o(980), 
the strong decay width is very large, T(/o — » 7T7t) = 
500 — 1000 MeV, according to various theoretical calcula- 
tions @. The small experimental width 40-100 MeV Q 
cannot be consistently explained by simple quark models. 

The strong-decay width suggests that /o(980) should 
not be an ordinary nonstrange gg-type meson. The 
Fermilab-E791 collaboration measured the decay D+ — > 
7r _ 7r + 7r + Q, which can proceed via intermediate states, 
for example, — > / (980)7r + with ss quarks in / (980). 
This experiment suggested a sizable strange-quark com- 
ponent in /o(980). The simplest configuration is a pure 
strange quarkonium ss for /o(980). In addition, since its 
mass is just below the K K threshold, it could be con- 
sidered as a KK molecule If two color-singlet states 
of K and K are not well separated, it corresponds to 
a tetraquark state, (uuss + ddss)/y/2, which was origi- 
nally suggested in the MIT bag model @. Recent QCD- 
sum-rule studies support this idea of a tetraquark state 
0- Furthermore, there are lattice-QCD studies that 



/o(980) corresponds to the tetraquark state because the 
scalar tetraquark mass is about 1.1 GeV @. In addition, 
/o(980) used to be considered as a glueball candidate; 
however, recent lattice QCD calculations rule out such 
a possibility because the mass of a ++ glueball is es- 
timated about 1700 MeV @. The situation of scalar 
mesons with masses in the 1 GeV region is summarized 
in Ref. [l[. All the possible /o(980) configurations are 
listed in Table Q] although the nonstrange-gg and glueball 
states seem to be unlikely according to the recent studies. 

In the following, the notation /o indicates the /o(980) 
meson and /o(600) is not discussed. There were propos- 
als to find the structure by a radiative decay into /o 
[1 , [ll], [l2| ■ Since it is an electric dipole decay, the width 
should reflect information on its size, namely its inter- 
nal structure [Io| . The experimental measurements of 
VEPP-2M [H[ and DA$NE [3 were reported for the 
decay <fi — ► /07. The data may suggest the tetraquark 
picture; however, there are still discussions on their in- 
terpretation [l2j]. Another possible experimental probe 
is the 77 — > 7r + 7r~ process in the /o mass region. The 
two-photon decay width of /o(980) was recently reported 
as 0.205^;^(stat)t£;"£(syst) keV by the Belle collabo- 
ration [15[]. Model calculations indicate 1.3 — 1.8 keV in 
the nonstrange gg picture; however, the measurement is 
consistent with the ss and KK -molecule configurations. 
There are also ideas to use elliptic flow and nuclear mod- 
ification ratios in heavy-ion reactions for finding exotic 
hadron structure fl6l ]. 

There are compelling theoretical and experimental ev- 
idences that the scalar meson /o(980) is not an ordinary 
nonstrange gg meson. However, a precise configuration 
is not determined yet, and a clear experimental evidence 
is awaited. It is the purpose of this paper to show that 
the internal structure of exotic hadrons should be de- 
termined from their fragmentation functions by noting 
differences in favored and disfavored functions. We in- 
vestigate Jo (980) as an example in this work. 

A fragmentation function is defined by a hadron- 



TABLE I: Possible /o(980) configurations and their features in fragmentation functions at small Q 2 . 



Type 


Configuration 


Second moments 


Peak positions 


Nonstrange qq 


(mm + dd) /V2 


M a <M U < Mg 


max max 
z u z s 


Strange qq 


ss 


M U <M S < Mg 


max , max 

z u z s 


Tetraquark (or KK) 


(uuss + ddss) /\/2 


M u ~M a <Mg 


max max 
z u z s 


Glueball 


99 


M u ~ M s < Mg 


max max 
z u ^ z s 



production cross section and the total hadronic cross sec- 
tion: F h (z,Q 2 ) = -L- d ^ e+e d ^ hx \ Here, the variable 
z is defined by the hadron energy Eh and the center- 
of-mass energy ^/s (= y%F) by z = Eh/{\fs/2). The 
fragmentation occurs from primary partons, so that it is 
expressed by the sum of their contributions: F h (z, Q 2 ) = 
J2 { Ci(z 7 a s )<S)D^(z, Q 2 ), where (g> indicates the convolu- 
tion integral, f(z)®g(z) = $\ dyf(y)g(z/y)/y, D^(z, Q 2 ) 
is the fragmentation function of the hadron h from a par- 
ton i (= u, d, s, <?), Ci(z,a s ) is a coefficient function, 
and a s is the running coupling constant. The favored 
fragmentation means a fragmentation from a quark or 
an antiquark which exists in a hadron as a constituent in 
a quark model, and the disfavored means a fragmentation 
from a sea quark. The favored and disfavored functions 
are assigned in the following discussions by considering 
the naive quark configurations in Table [TJ 

We first consider a possible ss configuration for fo. 
Then, the u- and <i-quark fragmentation functions are 
disfavored ones and the s-quark function is a favored 
one. For example, the favored fragmentation from s is 
possible if a gluon is radiated from s, and then it splits 
into a ss pair to form the /o meson as shown in Fig. [T] 
The notations 0(g 2 ) and 0(g 3 ) indicate the second and 
third orders of the coupling constant g. In the disfavored 
process from u, there are processes in the order of 0(g 3 ) 
without an 0(g 2 ) term, so that its probability is expected 
to be smaller than the favored one from s. It leads to the 
relation for the second moments of fragmentation func- 
tions: M u < M s , where M. t = JdzzD{°(z). The second 
moment Mj is the energy fraction for /o which is cre- 
ated from the parton i. In the same way, fragmentations 
occur from a gluon as shown in the figure. Since there 
are two processes in 0(g 2 ) with a soft gluon radiation, 
the second moment for the gluon is expected to be larger 
than the others. These considerations lead to the rela- 
tion M u < M s < M g in Table |U Such a naive estimation 
should be a crude one, but it has been shown to work for 
the moments of the pion, kaon, and proton [17| . so that 
it is also expected to be a reasonable guideline in other 
hadrons. 

Next, functional forms are discussed in the ss picture. 
More energy is transferred to fo from the initial s in the 
0(g 2 ) process than the one from the initial u due to an 
extra gluon emission in Fig. [T] It means that the frag- 
mentation function D[° (z) is distributed in the larger z 
region in comparison with D^ (z) because the /o energy 
is directly proportional to z. Namely, they should have 



different functional forms and their peak positions are 
different at small Q 2 (~1 GeV 2 ). We express this situ- 
ation as z™ ax < zf ax in Table [fl The form of the gluon 
fragmentation function may not be simply compared with 
the quark processes. 

In the same way, the second moments and functional 
forms are roughly estimated for the tetraquark picture. 
Since the fragmentations from u and s quarks are equally 
favored processes in this case, their moments and func- 
tions forms should be almost the same. The fragmen- 
tations into /o proceed by creating uu (or dd) and ss 
pairs as shown in Fig. [21 There are more fragmenta- 
tion processes from a gluon, so that the gluon moment 
is expected to be larger than the others. In this way, we 
obtain the relations, M u ~ M s < M g and zjf ax ~ z™ ax , 
in Table |TJ Since the flavor composition of fo is simply 
considered in the above discussions, this relation could 
be also applied to the KK case. However, the KK is 
a loose and extended bound state so that its production 
probability in the fragmentation is expected to be much 
smaller than that of the tetraquark state. 

Although the nonstrange-gg and glueball configura- 
tions seem to be unlikely according to recent theoreti- 
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FIG. 1: (Color online) Schematic diagrams for fo production 
in the ss configuration. 



Tetraquark picture for / (980) /; = ( U uss + ddss) 1 4l 
u, s (favored) 




+ six O (g 4 ) terms 
of gluon radiation 



0(g 4 ) 



FIG. 2: (Color online) Schematic diagrams for fo production 
in the tetraquark configuration. 
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cal investigations, we also estimated possible relations in 
Table [TJ Since the estimation method is essentially the 
same, derivations are not explained here. If /q were a 
nonstrange-gg meson, the relations M s < M M < M g and 
z™ ax > z™ ax are expected, whereas they are M u ~ M s < 
M g and z™ ax — z" lax if it were a glueball. 

The fragmentation functions are determined by a 
global analysis of hadron-production data in e + e~ an- 



nihilation [181 ]. There is recent progress on their anal- 
ysis. Uncertainties of the fragmentation functions are 
determined in Ref. [l7| , and it was shown that the gluon 
and light-quark functions have large uncertainties for the 
pion, kaon, and proton. Then, a global analysis with data 
in lepton scattering and proton-proton collisions was also 
reported This kind of global analysis is suitable for 
finding exotic hadrons by noting the typical features in 
the favored and disfavored functions. 

All the possible configurations for fo indicate that up- 
and down-quark compositions are the same; however, 
they are generally different from the strange-quark and 
other ones. Therefore, a natural and model-independent 
paramctrization is 



Dt(z,Q 2 Q ) = Dt(z,Ql) = D{\z,Ql) = D f °{z,Q 2 ) 

= Nl°z a t°(l-z)rt°, 
Df°(z,Q 2 ) = Di°{z,Ql) = AT/oX°(l - z fi\ 
D^{z,Ql)=N^z a i°{l-zfi\ (1) 
£>/o ( Z) ml) = Dt (z, m 2 c ) = N^z^ (1 - z)^° , 



D(°(z,mi) = Dl°(z,mi) = N'°z a i"(l - zf»\ 

where iVj, aj, and (3i are the parameters to be determined 
by a x 2 analysis of the data for e+ + e~ — > fo + X [2Cj |. 
The initial scale is taken Qq=1 GeV 2 , and the masses are 
m c =1.43 GeV and mb=4.3 GeV. The details of the anal- 
ysis method in the next-to-leading order are explained in 
Ref. [l7| • Uncertainties of the determined functions are 
estimated by the Hessian method [17J, which has been 
used also in the studies of various parton distribution 
functions EEL El: 



,2a _ n /o 



2 ^ _ KrfoW, 



[8D(°{z)Y = ^ 



' dD{°(z,Z) 



(2) 

Here, SDj" (z) is the uncertainty of the fragmentation 
function D{°(z), A% 2 value is taken so that the con- 
fidence level P becomes the one-a-error range (P = 
0.6826) by assuming the normal distribution in the multi- 
parameter space, Hij is the Hessian matrix, is a pa- 
rameter, and £ indicates the optimum parameter set. 

The number of fo data is very limited at this stage. In 
fact, available data are merely twenty three. This situ- 
ation makes the analysis difficult in obtaining the mini- 
mum \ 2 point. There are irrelevant parameters which do 
not affect the total % 2 . We decided to fix three parame- 
ters at /3 g = 1, ot u = 10, and a s = 10 because of the lack 
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FIG. 3: (Color online) Determined fragmentation functions 
of /o(980) by the global analysis. The functions zD(p , zD{° , 
and zD' g ° are shown at Q 2 =l GeV 2 , and the functions zD(° 
and zD(° are at Q 2 = m 2 and m 2 , respectively. 



of data. Then, the total number of parameters becomes 
twelve. The minimum x 2 is obtained x 2 /d.o.f.=0.907 in 
our analysis. 

The determined functions are shown in Fig. [3] It is 
interesting to find that the up- and strange-quark func- 
tions are distributed relatively at large z, and both func- 
tions have similar shapes, whereas the gluon, charm-, 
and bottom-quark functions are distributed at smaller 
z. It may indicate that both functions, D? and D{° , are 
equally favored ones, which implies that the up-quark 
(and down-quark) is one of main components of / as 
well as the strange-quark. Furthermore, they are peaked 
almost at the same points of z (z™ ax ~ z™ ax ), which 
may be also considered as an evidence for the tetraquark 
structure according to Table fl] However, if it is judged 
from their second moment ratio (M u /M s = 0.43), it looks 
like the ss configuration. 

This conflict is mainly caused by the inaccurate deter- 
mination of the fragmentation functions although it may 
be understood by admixture of the ss and tetraquark 
configurations. In Fig. 01 the uncertainties of zD^ , 
zD{°, and zD? are shown at Q 2 —l GeV 2 together 
with the functions themselves. We notice huge uncer- 
tainties which are an order of magnitude larger than 
the determined functions. If their moments are calcu- 
lated, they have large errors: M u = 0.0012 ± 0.0107, 
M s = 0.0027 ± 0.0183, and M g = 0.0090 ± 0.0046. From 
these results, the error of the moment ratio is estimated 
as M u /M s = 0.43 ± 6.73, which makes it impossible to 
discuss the effect of the order of 50%. In this way, we find 
the structure of fo cannot be determined by the current 
e + e~ data. 

It is the purpose of this work to point out that struc- 
ture of exotic hadrons should be determined by the frag- 
mentation functions. Accurate measurements of hadron- 
production cross sections can be used for determining 
their internal quark and gluon configurations as ex- 
plained in this paper by taking /o(980) as an example. 
We have shown that ss and tetraquark configurations, 
and also nonstrange-gg and glueball states, should be dis- 
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FIG. 4: (Color online) Fragmentation functions, zDl°, zD{° , 
and zDg° , and their uncertainties are shown at Q 2 =l GeV 2 . 
The uncertainties are shown by the shaded bands. 

tinguished by the second moments and functional forms 
of the favored and disfavored functions. Especially, the 
ratio of the u-quark moment to the s-quark one should 
be useful to judge the configuration. 

In order to determine the internal structure, the flavor 
separation is important especially because the difference 
between the up- and strange-quark functions is the key 
to find the structure of fo- First, charm- and bottom- 
quark tagged data should be provided for fo as they have 
been obtained for the pion, kaon, and proton. Then, 
the charm- and bottom-quark functions should be deter- 



mined accurately. Second, semi-inclusive /o-production 
data in lepton-proton scattering can be used for distin- 
guishing between up- and strange-quark fragmentations 
because the initial quark distributions arc different in the 
proton. These flavor separations will become possible 
by future experimental analyses. Our work is a starting 
point for exotic hadron search by suggesting the relations 
in the second moments and the functional forms and by 
indicating the current experimental situation as the un- 
certainty bands. 

The fragmentation functions of fo and their uncertain- 
ties have been determined by the global analysis of fo 
production data. At this stage, the e + e~ data are not 
precise enough; however, accurate experimental measure- 
ments could create a field of exotic hadrons which are be- 
yond the naive qq and qqq type ones. Currently, analyses 
are in progress by the Belle collaboration 23] to provide 
accurate fragmentation functions. They are especially 
important because the functions are measured at small 
Q 2 (^C Af§), so that scaling violation can be investigated 
to find the gluon functions [l7j . It is also important to 
have accurate measurements for ordinal mesons such as 
0(1020) and /2(1270) in order to establish the fo config- 
uration by comparing their favored and disfavored func- 
tions with the ones of fo- We could investigate other 
exotic hadrons in the same way by their fragmentation 
functions. 
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